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1. Introduction 

Neutron stars are dense, neutron-packed remnants of stars that blew apart in su- 
pernova explosions. Many neutron stars form radio pulsars, emitting radio waves 



With central densities way above the density of atomic nuclei, neutron stars contain 
matter in one of the densest forms found in the universe. Depending of the density 
reached in the cores of neutron stars, they may contain stable phases of exotic matter 
^ found nowhere else in space. This article gives a brief overview of the phases of ultra- 

I ' dense matter predicted to exist deep inside neutron stars and discusses the equation of 

O ' state (EoS) associated with such matter. 

> 
X 

^ i that appear from the Earth to pulse on and off like a lighthouse beacon as the 

star rotates at very high speeds. Neutron stars in x-ray binaries accrete material 
from a companion star and flare to life with a burst of x-rays. Measurements of 
radio pulsars and neutron stars in x-ray binaries comprise most of the neutron 
star observations. Improved data on isolated neutron stars (e.g. RX J1856. 5-3754, 
PSR 0205-1-6449) are now becoming available, and future investigations at gravi- 
tational wave observatories focus on neutron stars as major potential sources of 
gravitational waves (see Ref. 1 for a recent overview) . Depending on star mass and 
rotational frequency, the matter in the core regions of neutron stars may be com- 
pressed to densities that are up to an order of magnitude greater than the density of 
ordinary atomic nuclei. This extreme compression provides a high-pressure environ- 
ment in which numerous subatomic particle processes are believed to compete with 
each other^'^. The most spectacular ones stretch from the generation of hyperons 
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Fig. 1. Models for the EoS (pressure versus energy density) of neutron star matter^. The no- 
tation is as follows: RMF=relativistic mean-field model; DD-RBHF=density dependent rela- 
tivistic Brueckner-Hartree-Fock model; n=neutrons; p=protons; K=K~[u,s] meson condensate; 
quarks=u,d,s; H-mattcr=H-dibaryon condensate. 

and baryon resonances (E, A, S, A) to quark (u, d, s) deconfinement to the forma- 
tion of boson condensates (tt^, K^, H-mattcr)^''^'^'^'^'''. It has also been suggested 
(strange matter hypothesis) that strange quark matter may be more stable than 
ordinary atomic nuclei. In the latter event, neutron stars could in fact be made of 
absolutely stable strange quark matter rather than ordinary hadronic matter^'^'^". 
Another striking implication of the strange matter hypothesis is the possible exis- 
tence of a new class of white-dwarfs-like strange stars (strange dwarfs) The quark 
matter in neutron stars, strange stars, or strange dwarfs ought to be in a color su- 
perconducting state^^'^^'^^'^^. This fascinating possibility has renewed tremendous 
interest in the physics of neutron stars and the physics and astrophysics of (strange) 
quark matter^'-'^^'^^. This paper reviews the possible phases of ultra-dense nuclear 
matter expected to exist deep inside neutron stars and its associated equation of 

State2'3,4,5,6,7_ 

2. Neutron Star Masses 

In 1939, Tolman, Oppenheimer and Volkoff performed the first neutron star cal- 
culations, assuming that such objects are entirely made of a gas of relativis- 
tic neutrons"'^^'"'^^. The EoS of such a gas is extremely soft (i.e. very little addi- 
tional pressure is gained with increasing density), as can be seen from Fig. 1, 
and leads to a maximum neutron star mass of just 0.7 Mq. A relativistic neu- 
tron gas thus fails to accommodate neutron stars such as the Hulse- Taylor pul- 
sar (Af = 1.44 M0)^^, and also conflicts with the average neutron star mass of 
1.350 ± 0.004 Mq derived by Thorsett and Chakrabarty^^ from observations of 
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radio pulsar systems. More than that, recent observations indicate that neutron 
star masses may be as high as ~ 2 Mq. Examples of such very heavy neu- 
tron stars arc Aif,io75i+i807 = 2.1 ±0.2 Mq^", Aif4u 1636+536 = 2.0 ± 0.1 Mq^i, 
Mvciax-i = 1.86 ± 0.16 M022, Mcygx-2 = 1.78 ± 0.23 Aif023.24^ Large masses have 
also been reported for the high-mass x-ray binary 4U 1700-37 and the compact ob- 
ject in the low-mass x-ray binary 2S0921-630, M4U 1700-37 = 2.44 ±0.27 Mg^^ and 
-^^230921-630 = 2.0 — 4:.3Mq'^^ . rcspcctivcly. The latter two objects may be cither 
massive neutron stars or low-mass black holes with masses slightly higher than the 
maximum possible neutron star mass of ~ 3Mq. This value follows from a general, 
theoretical estimate of the maximal possible mass of a stable neutron star as per- 
formed by Rhoadcs and Ruffini^"^ on the basis that (1) Einstein's theory of general 
relativity is the correct theory of gravity, (2) the EoS satisfies both the microscopic 
stability condition dP/de > and the causality condition dP/de < c^, and (3) 
that the EoS below some matching density is known. From these assumptions, it 
follows that the maximum mass of the equilibrium configuration of a neutron star 
cannot be larger than 3.2 Mq. This value increases to about 5Mq if one abandons 
the causality constraint dP/de < c^,^^'^^ since it allows the EoS to behave stiffer at 
asymptotically high nuclear densities. If either one of the two objects 4U 1700-37 
or 2S092 1-630 were a black hole, it would confirms the prediction of the existence 
of low-mass black holes^°. Conversely, if these objects were massive neutron stars, 
their high masses would severely constrain the EoS of nuclear matter. 

3. Composition of Ultra-Dense Neutron Star Matter 

Models for the EoS of neutron star matter are being computed in different 
theoretical frameworks. The most popular ones are the semi-classical Thomas- 
Fermi theory^-'^'^^^ Schroedinger-based treatments (e.g. variational approach, Monte 
Carlo techniques, hole line expansion (Brueckner theory), coupled cluster method. 
Green function method)^'^^'^^'^^, or relativistic field-theoretical treatments (rela- 
tivistic mean field (RMF), Hartree-Fock (RHF), standard Brueckner-Hartree-Fock 
(RBHF), density dependent RBHF (DD-RBHF)36^37,38,39404i^ the Nambu-Jona- 
Lasinio (NJL) model'*^''*^, and the chiral SU(3) quark mean field model^'*). An 
overview of EoSs computed for several of these methods is given in Fig. 1. Mass- 
radius relationships of models of neutron stars based on these EoSs are shown in 
Fig. 2. Any acceptable nuclear many-body calculation must reproduce the bulk 
properties of nuclear matter at saturation density and of finite nuclei. The nuclear 
matter properties are the binding energy, E/A = —16.0 MeV, effective nucleon 
mass, mj^ = 0.79 mN, incompressibility, K ~ 240 MeV, and the symmetry energy, 
ttg = 32.5 MeV, at a saturation density of no = 0.16 fm~^. 

3.1. Hyperons 

Only in the simplest conception, a neutron star is constituted from only neutrons. 
As already discussed in Sect. 2, this model fails by far to accommodate observed 
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Fig. 2. Mass-radius relationship of neutron stars and strange stars. The strange stars may be en- 
veloped in a crust of ordinary nuclear material whose density is below neutron drip density*'*^'*^ . 
The labels are explained in Fig. 1. (Figure from Ref. 6.) 



neutron star masses. At a more accurate representation, neutron stars will contain 
neutrons, n, and a small number of protons, p, whose charge is balanced by elec- 
trons, e^, according to the reaction n p + e^ +Ve and its inverse. At the densities 
that exist in the interiors of neutron stars, the neutron chemical potential, eas- 
ily exceeds the mass of the A so that neutrons would be replaced with A hyperons. 
From the threshold relation /i" — fi^ it follows that this would happen for neutron 
Fermi momenta greater than kp^ 3fm^^.^ Such Fermi momenta correspond to 
densities of just ~ 2no. Hence, in addition to nuclcons and electrons, neutron stars 
may be expected to contain considerable populations of strangeness-carrying A hy- 
perons, possibly accompanied by smaller populations of the charged states of the 
S and S hyperons*^^. The total hyperon population may be as large as 20% ^'^ . 

3.2. Meson condensation 

The condensation of negatively charged mesons in neutron star matter is favored 
because such mesons would replace electrons with very high Fermi momenta. Early 
estimates predicted the onset of a negatively charged pion condensate at around 
2no (see, for instance, Ref. 48). However, these estimates are very sensitive to the 
strength of the effective nucleon particle- hole repulsion in the isospin T = 1, spin 
S = 1 channel, described by the Landau Fermi- liquid parameter g', which tends 
to suppress the condensation mechanism. Measurements in nuclei tend to indicate 
that the repulsion is too strong to permit condensation in nuclear matter^^''"^". In 
the mid 1980s, it was discovered that the in-medium properties of K~ [us] mesons 
may be such that this meson rather than the 7r~ meson may condense in neutron 
star matter^^'^^'^^. 
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The condensation is initiated by the schematic reaction e~ — > K~ + Ug. If this 
reaction becomes possible in neutron star matter, it is energetically advantageous 
to replace the fermionic electrons with the bosonic K~ mesons. Whether or not 

this happens depends on the behavior of the K~ mass, 'Ti^_ , in neutron star mat- 
ter. Experiments which shed light on the properties of the K~ in nuclear matter 
have been performed with the Kaon Spectrometer (KaoS) and the FOPI detector 
at the heavy-ion synchrotron SIS at GSi54,55,56,57,58 analysis of the early K~ 
kinetic energy spectra extracted from Ni+Ni collisions'"'^ showed that the attraction 
from nuclear matter would bring the mass down to 2± 200 MeV at den- 

sities ~ 3no. For neutron-rich matter, the relation m'^- /rrix- — 1 — 0.2n/no was 
established^^'®"'®^, with niK = 495 MeV the K~ vacuum mass. Values of around 
~ 200 MeV may be reached by the electron chemical potential, n'^, in neutron 
star matter^'^^ so that the threshold condition for the onset of K~ condensation, 
1^^ — might be fulfilled for sufficiently dense neutron stars, provided other 
negatively charged particles A~, d and s quarks) are not populated first and 
prevent the electron chemical potential from increasing monotonically with density. 

We also note that K~ condensation allows the conversion reaction n p + K~. 
By this conversion the nucleons in the cores of neutron stars can become half 
neutrons and half protons, which lowers the energy per baryon of the matter^^. 
The relative isospin symmetric composition achieved in this way resembles the one 
of atomic nuclei, which are made up of roughly equal numbers of neutrons and 
protons. Neutron stars are therefore referred to, in this picture, as nucleon stars. 
The maximum mass of such stars has been calculated to be around 1.5 Mq^^. 
Consequently, the collapsing core of a supernova, e.g. 1987A, if heavier than this 
value, should go into a black hole rather than forming a neutron star, as pointed 
out by Brown et al.'^"'^^'^" This would imply the existence of a large number of 
low- mass black holes in our galaxy'^°. Thielemann and Hashimoto^^ deduced from 
the total amount of ejected ^^Ni in supernova 1987A a neutron star mass range 
of 1.43 — 1.52 Mq. If the maximum neutron star mass should indeed be in this 
range (~ 1.5 Mq), it would pose a problem to the possible discovery of very heavy 
neutron stars of masses around 2 Mq (Sect. 2). In closing, we mention that meson 
condensates lead to neutrino luminosities which are considerably enhanced over 
those of normal neutron star matter. This would speed up neutron star cooling 
considerably^^'^^. 

3.3. H-matter and exotic baryons 

A novel particle that could be of relevance for the composition of neutron star mat- 
ter is the H-dibaryon (H=(['ud][ds][su])), a doubly strange six-quark composite with 
spin and isospin zero, and baryon number two^^. Since its first prediction in 1977, 
the H-dibaryon has been the subject of many theoretical and experimental studies 
as a possible candidate for a strongly bound exotic state. In neutron star matter, 
which may contain a significant fraction of A hyperons, the A's could combine to 
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form H-dibaryons, which could give way to the formation of H-dibaryon matter at 
densities somewhere above ~ 4no^'^'^^'^^. If formed in neutron stars, however, H- 
matter appears to unstable against compression which could trigger the conversion 
of neutron stars into hypothetical strange stars^^'™''^-'^. 

Another particle, referred to as exotic baryon, of potential relevance for neutron 
stars, could be the pentaquark, ©+([ud]^s), with a predicted mass of 1540 MeV. 
The pentaquark, which carries baryon number one, is a hypothetical subatomic 
particle consisting of a group of four quarks and one anti-quark (compared to three 
quarks in normal baryons and two in mesons), bound by the strong color-spin 
correlation force (attraction between quarks in the color 3c channel) that drives 
color superconductivity^^'^''. The pentaquark decays according to O+(1540) 
K~^[su\ + n[udd] and thus has the same quantum numbers as the X+n. The asso- 
ciated reaction in chemically equilibrated matter would imply /x®^ = ij,^* + n^. 



3.4. Quark deconfinement 

It has been suggested already many decades ago'''^'^^''''^'^^''^®''^^'^'' that the nucleons 
may melt under the enormous pressure that exists in the cores of neutron stars, 
creating a new state of matter know as quark matter. From simple geometrical 
considerations it follows that for a characteristic nucleon radius of r^v ~ 1 fm, 
nucleons may begin to touch each other in nuclear matter at densities around 
(47rr^/3)~^ ~ 0.24 fm^'^ = 1.5 no, which is less than twice the density of nuclear 
matter. This figure increases to ~ lino for a nucleon radius of rjv = 0.5 fm. One 
may thus speculate that the hadrons of neutron star matter begin to dissolve at 
densities somewhere between around 2 — 10 no, giving way to unconfined quarks. 
Depending on rotational frequency and neutron star mass, densities greater than 
two to three times no are easily reached in the cores of neutron stars so that the 
neutrons and protons in the cores of neutron stars may indeed bo broken up into 
their quarks constituents^'^'^'^^. More than that, since the mass of the strange 
quark is only nig ~ 150 MeV, high-energetic up and down quarks will readily 
transform to strange quarks at about the same density at which up and down 
quark deconfinement sets in. Thus, if quark matter exists in the cores of neutron 
stars, it should be made of the three lightest quark flavors. A possible astrophysical 
signal of quark deconfinement in the cores of neutron stars was suggested in Ref. 82. 
The remaining three quark flavors (charm, top, bottom) are way to massive to be 
created in neutron stars. For instance, the creation of charm quark requires a density 
greater than lO^^g/cm^, which is ~ 10^ times greater than the density reached in 
neutron stars. A stability analysis of stars with a charm quark population reveals 
that such objects are unstable against radial oscillations and, thus, can not exist 
stably in the universe^'®. The same is true for ultra-compact stars with unconfined 
populations of top and bottom quarks, since the pulsation eigen-equations are of 
Sturm-Liouville type. 
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Fig. 3. Radii of quark bags For masses less than 10^ GeV the electrons (gray dots) are outside 
the quark bags (indicated by thick solid circles) and the core+electron system has a size of ~ 
10^ fm. For masses between 10^ and 10^® GeV the electrons are partially inside the core. For 
masses greater than 10^^ GeV all electrons are inside the core, except for a small electron fraction 
which forms an electric dipole layer on the surface of the quark bag (strange star). 

4. Strange Quark Matter 

It is most intriguing that for strange quark matter made of more than a few hundred 

up, down, and strange quarks, the energy of strange quark matter may be well 
below the energy of mielear matter, E/A = 930 MeV, which gives rise to new and 
novel classes of strange matter objects (see Fig. 3), ranging from strangelets at 
the low baryon-number end to strange stars at the high baryon number end. A 
simple estimate indicates that for strange quark matter E/A = ABw^ /n^, so that 
bag constants of B = 57 MeV/fm^ (i.e. B^^"^ = 145 MeV) and 5 = 85 MeV/fm^ 
(B^/^ = 160 MeV) would place the energy per baryon of such matter at E/A = 
829 MeV and 915 MeV, respectively, which correspond obviously to strange quark 
matter which is absolutely bound with respect to nuclear matter*"'' 

4.1. Nuclear crust on strange stars 

Strange quark matter is expected to be a color superconductor which, at extremely 
high densities, should be in the CFL phase^^'^^. This phase is rigorously electrically 
neutral with no electrons required*^. For sufficiently large strange quark masses, 
however, the low density regime of strange quark matter is rather expected to form 
other condensation patterns (e.g. 2SC, CFL-if^, CFL-K+, CFL-7r° ") in which 
electrons are present^^'^^. The presence of electrons causes the formation of an 
electric dipole layer on the surface of strange matter, with huge electric fields on 
the order of 10^^ V/cm, which enables strange quark matter stars to be enveloped 
in nuclear crusts made of ordinary atomic matter*'^'*^'*^.*^ The maximal possible 

^Depending on the surface tension of blobs of strange matter and screening effects, a heterogeneous 
crust comprised of blobs of strange quark matter embedded in an uniform electron background 
may exist in the surface region of strange stars** . This heterogeneous strange star surface would 
have a negligible electric field which would make the existence of an ordinary nuclear crust, which 
requires a very strong electric field, impossible. 



5, 2008 11:37 WSPC/INSTRUCTION FILE ishipTWeber 



8 F. Weber et al. 




Density (g cm ) 



Fig. 4. Illustration of the EoS of strange stars with nuclear crusts (from Ref. 86) . 



density at the base of the cnist (inner crust density) is determined by neutron drip, 
which occurs at about 4 x 10^^ g/cni' or somewhat below^^. The EoS of such a 
system is shown in Fig. 4. Sequences of compact strange stars with and without 
(bare) nuclear crusts are shown in Fig. 2. Since the nuclear crust is gravitationally 
bound to the quark matter core, the mass-radius relationship of strange stars with 
crusts resembles the one of neutron stars and even that of white dwarfs^^. Bare 
strange stars obey M (x E? because the mass density of quark matter is almost 
constant inside strange stars. 



4.2. Strange dwarfs 

For many years only rather vague tests of the theoretical mass-radius relationship of 
white dwarfs were possible. Recently the quality and quantity of observational data 
on the mass-radius relation of white dwarfs has been reanalyzed and profoundly 
improved by the availability of Hipparcos parallax measurements of several white 
dwarfs^^. In that work Hipparcos parallaxes were used to deduce luminosity radii 
for 10 white dwarfs in visual binaries of common proper-motion systems as well as 
11 field white dwarfs. Complementary HST observations have been made to better 
determine the spectroscopy for Procyon B^" and pulsation of G226-29^^. Procyon B 
at first appeared as a rather compact star which, however, was later confirmed to 
lie on the normal mass-radius relation of white dwarfs. Stars like Sirius B and 
40 Erin B, fall nicely on the expected mass-radius relation too. Several other stars 
of this sample (e.g. GD 140, G156-64, EG 21, EG 50, G181-B5B, GD 279, WD2007- 
303, G238-44) however appear to be unusually compact and thus could be strange 
dwarf candidates The situation is graphically summarized in Fig. 5. 
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Fig. 5. Comparison of the theoretical mass-radius relationships of strange dwarfs (solid curves) 
and normal white dwarfs®®. Radius and mass are in units of Rq and Mq, respectively. 



4.3. Surface properties of strange matter 

The electrons surrounding strange quark matter are held to quark matter elec- 
trostatically. Since neither component, electrons and quark matter, is held in place 
gravitationally, the Eddington limit to the luminosity that a static surface may emit 
does not apply, and thus the object may have photon luminosities much greater than 
10^* erg/s. It was shown by Usov^^ that this value may be exceeded by many orders 
of magnitude by the luminosity of e+e" pairs produced by the Coulomb barrier at 
the surface of a hot strange star. For a surface temperature of ^ 10^^ K, the Imni- 
nosity in the outflowing pair plasma was calculated to be as high as ^ 3 x 10°^ erg/s. 
Such an effect may be a good observational signature of bare strange stars^'^'^'^'^^'^^. 
If the strange star is enveloped by a nuclear crust however, which is gravitationally 
bound to the strange star, the surface made up of ordinary atomic matter would 
be subject to the Eddington limit. Hence the photon emissivity of such a strange 
star would be the same as for an ordinary neutron star. If quark matter at the 
stellar surface is in the CFL phase the process of e+e~ pair creation at the stellar 
quark matter surface may be turned off, since cold CFL quark matter is electrically 
neutral so that no electrons are required and none are admitted inside CFL quark 
matter*^. This may be different for the early stages of a hot CFL quark star^''. 

5. Proto-Neutron Star Matter 

Here we take a brief look at the composition of proto-neutron star matter. The 
composition is determined by the requirements of charge neutrality and equilib- 
rium under the weak processes, Bi ^ B2 + I + Pi and B2 + 1 ^ Bi + i^i, where Bi 
and B2 are baryons, and I is a lepton, either an electron or a muon. For standard 
neutron star matter, where the neutrinos have left the system, these two require- 
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ments imply that Q — J2i 'Zi^-Si + J2i=e ^ ^I'^i ~ ^ (electric charge neutrahty) and 
^■Bi = bifin — QifJ-i (chemical equilibrium), where qi/i denotes the electric charge 
density of a given particle, and (ni) is the baryon (lepton) number density. The 
subscript i runs over all the baryons considered. The symbol fiBi refers to the chem- 
ical potential of baryon i, bi is the particle's baryon number, and Qi is its charge. 
The chemical potential of the neutron is denoted by /x„. When the neutrinos are 
trapped, as it is the case for proto-neutron star matter, the chemical equilibrium 
condition is altered to /x^^ = bi^n — <li(,l^i — ^J^Vl) and — ^jl^^ = IJi-n — fJ-i^,^, where 
fi^i is the chemical potential of the neutrino i>i. In proto-neutron star matter, the 
electron lepton number — (rie -I- )/nB is initially fixed at a value of around 
Yl^ = Ye + Yy^ ~ 0.3 — 0.4 as suggested by gravitational collapse calculations of 
massive stars. Also, because no muons are present when neutrinos are trapped, the 
constraint Yl^ = Y^j_ + Y^^u = can be imposed. Figures 6 and 7 show sample 
compositions of proto-neutron star matter and standard neutron star matter (no 
neutrinos) computed for the RMF approximation. The presence of the A particle 
in (proto) neutron star matter at finite temperature is striking. This particle is 
generally absent in cold neutron star matter treated in RMF^''^'^^. 

6. Electrically Charged Relativistic Star 

As described in Sect. 4.1, the electric field existing on the surface of a strange 
quark star could be as high as 10^^ V/cm^^. The energy density associated with 
such a ultra-high field begins to have an impact one the geometry of space-time, 
as determined by Einstein's field equation^^"'^'^^'^"^. In what follows, we investigate 
the effects of such ultra-high electric fields on the properties of relativistic polytropic 
stars. The consequences for strange stars are explored in Ref. 103. Since we want 
to maintain spherical symmetry of the star, the natural choice for the metric is 

ds^ = e^'^^^c^dt' - e^^'-^dr^ - r\de^ + sin^ Odcf") . (1) 
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Fig. 8. Star's mass existing as baryonic mass, 
mt, and electric mass-energy, me, for / = 
0.0008. 



Fig. 9. Same as Fig. 8, but for / = 0.001 



The energy-momentum tensor of the star consists of two terms, the standard term 
which describes the star's matter as a perfect fluid and the electromagnetic term, 



T,'' = {p + pc^)u^ u^+pS^^' + l- (^F^'^F^i + i- 5," Fki F'^') , 



(2) 



where the components F'^'^ satisfy the covariant Maxwell equations [{—9Y^'^F'^^],pi 
= 47rJ'^(— 5)^/^, with J" the four current. For the radial component one obtains 

(3) 

(4) 
(5) 



F°\r) = E{r) = e-(''+^)/2 r'^ / A^fe^'^+^'^/Hr' . 

Jo 

From last equation we can define the charge of the system as 

Q{r)= r A-Kfr'^e^^+^y^dr'. 
Jo 

The electric field is then given by 

E{r) = 6-^"+''^^^ r-^ Q{r) . 
With the aid of these relations the energy-momentum tensor takes the form 



rp IJ, 



V 


























(6) 



Substituting this relation into Einstein's field equation, G^'^ = (87rG/c^)T^'', leads 
to 



l_ _ SttG / _ Q^(r) 
' r dr J \ Snr'^ 



1 IdX 

n + 

r 



1 dv 
r dr 



1 



1 



(7) 

(8) 
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Table 1. Hcsulls 



/ M (Mq) Mi (Mq) Me (Mq) Ratio Charge (xlO" C) 

0.0 1.428 1.428 0% 

0.0001 1.439 1.436 0.00250 0.17% 260.95 

0.0008 2.548 2.166 0.38227 15% 3488.01 

0.001 4.156 2.905 1.250 30.07% 6727.17 



The solution for the metric function A(r) is given by 



The first two terms on the right-hand-side of Eq. (9) describe electrically neutral 
stars, while the third term originates from the net electric charge distribution inside 
the star. Prom Eqs. (7) (9) one finds that the mass contained in a spherical shell 
of radius, m(r), is given by 

dmjr) _ Airr^ ^ ^ Qjr) dQ{r) 
dr (?r dr 

The first term on the right-hand-side is the standard result for the gravitational 
mass of electrically uncharged stars, while the second term accounts for the mass 
change that originates from the electric field. Finally, the Tolman-Oppenheimer- 
Vokoff (TOV) equation generalized to electrically charged stars follows as 



dp 



dr ^2^2 A _ 2Gm(r) , GQ^r)\ ^ 4nr^ dr 



(p+^) + '^^- (11) 



To solve the TOV equation, we need to specify the charge distribution inside the 
star. Here, we will follow the approach of Ref. 100, 104 and assume that the charge 
distribution is proportional to the energy density, = / x e, where / is a 

constant which essentially controls the amount of net electric charge carried by the 
star. Moreover, we choose as initial and boundary conditions of the problem the 
following values, e(0) = 1550 MeV/fm^, Q(0) = m(0) = A(0) = !/(0) = 0, and 
p{R) = 0. In this study we focus on the difference between the mass-energy coming 
from baryonic matter, here represented by the polytropic EoS, and from the electric 
field. To this aim we rewrite Eq. (10) as (6=baryonic matter, e=electric field) 

dm{r) dnib drue ^ drrib Anr"^ driie Q(r) dQ{r) 
dr dr dr ' rfr ' dr (?r dr 

Table 6 summarizes the propcirties of several cilcictrically charged sample stars. Fig- 
ures 8 and 9 show, for two selected values of /, how much of the star's total (grav- 
itational) mass exists as baryonic mass (mb) and how much mass-energy (rrie) is 
associated with the electric field energy. For a strongly charged star, / = 0.001, we 
find that 30% of the star's mass (1.25 Mq) is actually electrostatic energy. 
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